On Subsystem Codes Beating the 
Hamming or Singleton Bound 

By Andreas Klappenecker and Pradeep Kiran Sarvepalli 
Texas A&M University, College Station, TX 77843, USA 

Subsystem codes are a generalization of noiseless subsystems, decoherence free sub- 
spaces, and quantum error-correcting codes. We prove a Singleton bound for ¥ q - 
linear subsystem codes. It follows that no subsystem code over a prime field can 
beat the Singleton bound. On the other hand, we show the remarkable fact that 
there exist impure subsystem codes beating the Hamming bound. A number of 
open problems concern the comparison in performance of stabilizer and subsystem 
codes. One of the open problems suggested by Poulin's work asks whether a subsys- 
tem code can use fewer syndrome measurements than an optimal MDS stabilizer 
code while encoding the same number of qudits and having the same distance. We 
prove that linear subsystem codes cannot offer such an improvement under complete 
fT) , decoding. 
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O ' 1. Introduction 



Subsystem codes (sometimes also referred to as operator quantum error-correcting 
codes) have emerged as an important new discovery in the area of quantum error 
correcting codes, unifying the classes of stabilizer codes, decoherence free subspaces 
and noiseless subsystems (Bacon 2006; Knill 2006; Kribs et al. 2005, 2006; Kribs 
2006; Poulin 2005). From a practical perspective their importance lies in the fact 
that they seem to offer better error recovery schemes than existing quantum codes. 
Therefore, it is crucial to know under what circumstances these gains can be at- 
tained and how to achieve them. 

Recall that a quantum code Q is a subspace in a finite dimensional Hilbert space, 
H = C q .A subsystem code is a quantum code which can be further resolved into a 
tensor product i.e., Q — A® B . Information is stored in system A, while system B, 
referred to as the gauge subsystem, provides some additional redundancy. By qudit 
we refer to a quantum bit with q levels. We denote the parameters of a subsystem 
code by [[n, k, r, d\] q , indicating that it is a g-ary code with length n, encodes k 
qudits into the subsystem A, and contains r gauge qudits and has distance d. 

Our goals in this paper are twofold. After reviewing the necessary background 
on subsystem codes, we generalize the quantum Singleton bound to F g -linear sub- 
system codes. It follows that no Clifford subsystem code over a prime field can beat 
the Singleton bound. We use these results to show that if there exists an MDS 
stabilizer code, then no linear subsystem code can outperform it in the sense of 
requiring fewer syndrome measurements for error correction. 
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Bacon & Casaccino (2006) obtain a subsystem code from two classical codes. We 
show that this method is a special case of the Euclidean construction for subsystem 
codes proposed in Aly et al. (2006) and give a coding theoretic analysis of these 
codes. 

Since the early works on quantum error-correcting codes, it has been suspected 
that impure codes should somehow perform better than the pure codes. In particu- 
lar, it was often conjectured that there might exist impure quantum error-correcting 
codes beating the quantum Hamming bound, but a proof remained elusive. Aly et 
al. (2006) proved a Hamming bound for pure subsystem codes. We show here that 
there exist impure subsystem codes beating the Hamming bound. 

2. Background 

Let F q be a finite field with q elements and characteristic p. Let C C F™ be an ¥ q - 
linear classical code denoted by [n, k, d] q , where k = dimp^ C and d is the minimum 
distance of C. We define wt(C) = min{wt(c) / c e C} = d, where wt(c) is 
the Hamming weight of c. Sometimes an alternative notation (n, K, d) q is also used 
where K = \C\. If C is an F p -linear subspace over F g , then we say it is an additive 
code. 

If x, y € F™, then their Euclidean inner product is defined as x- y — J2i x iV%- The 
Euclidean dual of a code C C is defined as C 1 - = {y e ¥ q | x ■ y = for all x e 
C}. We say that a code C is self-orthogonal with respect to the Euclidean inner 
product ifCCC 1 . 

We use the notation (x\y) = (xi, . . . , x n \yi, . . . , y n ) to denote concatenation of 
x, y S F™. Let u = (a\b) and v — (a'\b') be in F^™. We define the symplectic weight of 
u as swt(u) = {(dj, bi) ^ (0, 0) | 1 < i < n} and the symplectic weight of a code C C 
¥ 2 q n as swt(C) = min{swt(c) ^ c e C}. For codes over F^™ another inner product 
plays a more important role in the context of quantum codes. The trace-symplectic 
product between u,v is defined as (u\v) t = {(a\b)\(a'\b')) t = tr q / p (a' ■ b — a ■ b'). 
The trace-symplectic dual of C C ¥ 2 q n is defined as C u = {x e ¥ 2 q n \ (x\y) t = 
0, for all y G C}. If C C C^*, we say that it is self-orthogonal with respect to the 
trace-symplectic inner product. 

(a) Subsystem codes from classical codes 

We now briefly review the background on subsystem codes. First we give a group 
theoretic description and then give an alternate description in terms of classical 
codes. Further details can be found in Klappenecker & Sarvepalli (2006); Aly et al. 
(2006). 

Let q be the power of a prime p and ¥ q a finite field with q elements. Let 
B = {|x) | x S F g } denote an orthonormal basis for C q . Let X(a) and Zip) be 
unitary operators on C q whose action on any element |x) in B is defined as 

X(a) |x) = \x + a) and Zip) |x) = cu tr "/" (bx) \x) , 

where uj = eJ 27T / p is a primitive p th root of unity. These operators are a g-ary gen- 
eralization of the well-known Pauli matrices X and Z. Their action on an arbitrary 
element in C 9 is obtained by invoking linearity Let TL — C q ® • • • ® C q = C 9 " and 
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£ be the error group on H, denned as the tensor product of n such error operators 
i.e., 

£ = {lo c E 1 ®---®E n \Ei = X(ai)Z{bi)\ a h h £¥ q ;c£ ¥ p }. 

The weight of an error E = lo c E\ ® E2 <S> ■ ■ ■ <S> E n in £ is defined as the number of Ei 
which are not equal to identity and it is denoted by wt(E). We can also associate 
to£a vector E — (ai, . . . , a n \bi, . . . , b n ) £ ¥ 2n . We define the symplectic weight of 
E as 

swt(S) = \{(a,i, bi) ^ (0, 0) I 1 < i < n}\ = wt(E). 

Every nontrivial normal subgroup N in £ defines a subsystem code Q. Let 
Cg(N) be the centralizer of N in £ and Z(N) the center of N. As a subspace the 
subsystem code Q defined by N is precisely the same as the stabilizer code defined 
by Z(N). By Theorem 4 in Klappenecker & Sarvepalli (2006), Q can be decomposed 
as A ® B where dimS = \G : Z{N)\ 1 / 2 and 

dirndl = \Z{£) n G\\£ : Z^ 1 ' 2 ^ : Z{N)\^ 2 /\N\. 

Since information is stored only on subsystem A, we need only concern errors that 
affect A. An error E in £ is detectable by subsystem A if and only if E is contained 
in the set £ — {NCg{N) — N). The distance of the code is defined as 

d = min{wt(£0 \ I^E£ NC £ (N) - N} = wt(NC £ (N) ~ N). 

If NCs(N) = N, then we define the distance of the code to be wt(N). A distance 
d subsystem code with dirndl = K, dim_B = R is often denoted as ((n, K, R,d)) q 
or [[n, k, r, d]] q if K = q k and R = q r . We say that N is the gauge group of Q and 
Z(N) its stabilizer. The gauge group acts trivially on A. 

In Klappenecker & Sarvepalli (2006) we showed that subsystem codes, much 
like the stabilizer codes, are related to the classical codes over F^" or F™ 2 , but with 
one important difference. We no longer need the associated classical codes to be 
self-orthogonal, thereby extending the class of quantum codes. The gauge group 
N can be mapped to a classical code C over ¥ 2n and Ce{N) can be mapped to 
the trace-symplectic dual of C . The following theorem (Klappenecker & Sarvepalli 
2006) shows how subsystem codes are related to classical codes . 

Theorem 2.1. Let C be a classical additive subcode of F^™ such that C 7^ {0} and 
let D denote its subcode D = Cf)C *. If x = \C\ and y — \D\, then there exists an 
operator quantum error correcting code C — A® B such that 
1) dim A = q n /{xy) 1 ' 2 , 
ii) dimB = (x/y) 1 / 2 . 

The minimum distance of subsystem A is given by 

(a) d = swt((C + C u ) - C) = swt(D- L * - C) if ^ C; 

(b) d = swt(D^) if = C. 

Thus, the subsystem A can detect all errors in £ of weight less than d, and can 
correct all errors in E of weight < [(d — l)/2j . 

We call codes constructed using theorem 2.1 as Clifford subsystem codes. Ar- 
guably, these codes cover the most important subsystem codes, including the re- 
cently proposed Bacon-Shor codes. In this paper, henceforth by a subsystem code 
we will mean a Clifford subsystem code. 
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A further simplification of the above construction is possible which takes any 
pair of classical codes to give a subsystem code. We will just recall the result here 
and study its application in the next section. 

Corollary 2.2 (Euclidean Construction). Let X, C F™ ; be [n,ki] q linear codes 
where i E {1,2}. Then there exists an [[n,k,r,d]] q Clifford subsystem code with 

• k = n - (fci + k 2 + k')/2, 

• r = (k 1 +k 2 - k')/2, and 

• d = min{wt((X 1 J - n X 2 )^ \ Xr), wt((X 2 1 - n Xi) x \ X 2 )}, 
w/iere fe' = dim Fg (Xi n X^) x (Xf n X 2 ). 

The result follows from Theorem 2.1 by defining C — X\ x X 2 ; it follows that 
C- 1 - = X^ x Xf andD^Cn C" 1 * = (Xi nl^) x (X 2 nXf), and the parameters 
are easily obtained from these definitions, see Aly et al. (2006) for a detailed proof. 

(b) Pure and impure subsystem codes 

We can extend the notion of purity to subsystem codes also in a straightforward 
manner. Let N be the gauge group of a subsystem code Q with distance d = 
wt(Cf(Z(iV)) — N). We say that Q is pure to d' if there is no error of weight less 
than d! in N. The code is said to be exactly pure to d' if wt(iV) is d' and it is said 
to pure if d' > d . The code is said to be impure if it is exactly pure to d' < d. This 
refinement to the notion of purity was made in recognition of certain subtleties that 
had to addressed when constructing other subsystem codes from existing subsystem 
codes, see Aly et al. (2006) for details. 

In coding theoretic terms this can be translated as follows. Let C be an additive 
subcode of V 2 q n and D = C n C u . By theorem 2.1, we can obtain an ((n, K, R,d)) q 
subsystem code Q from C that has minimum distance d — swt(D ±t — C). If d' < 
swt(C), then we say that the associated operator quantum error correcting code is 
pure to d! . 

Extending these ideas of purity to subsystem codes is useful because it facilitates 
the analysis of the parameters of the subsystem codes, as will become clear when 
we derive bounds in the next section. If the codes are pure, then it will be very easy 
to see that the subsystem code with the parameters [[n, k, r, d]] q satisfies k + r < 
n — 2d + 2. This is because then the subsystem code can also be viewed as an 
[[n, k + r, d]] q stabilizer code, see theorem 11 in Aly et al. (2006) for further details. 

3. Singleton upper bound for F^-linear subsystem codes 

(a) An upper bound for subsystem codes 

We prove that the F g -linear subsystem codes with the parameters [[n, k, r, d]] q 
satisfy a quantum Singleton like bound viz., k + r < n — 2d + 2. It will be seen 
that this reduces to the quantum Singleton bound if r = 0. More interestingly, this 
reveals that there is a trade off in the size of subsystem A and the gauge subsystem. 
One pays a price for the gains in error recovery. The cost is the reduction in the 
information to be stored. 

Our proof for this result is quite straightforward, though the intermediate details 
are a little involved. First we show that a linear [[n, k, r > 0, d]] q subsystem code 
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that is exactly pure to 1 can be punctured to an [[n — 1, k, r — 1, d]] q code which 
retains the relationship between n, k, r, d. If d = 2 by repeated puncturing we either 
arrive at a pure code or a stabilizer code, both of which have upper bounds. For 
d > 2, two cases can arise, if the code is exactly pure to 1, we simply puncture 
it to get a smaller code as in d = 2 case. Otherwise, we puncture it to get an 
[[n — 1, k, r + 1, d — l]] q code. By repeatedly shortening we either get a stabilizer 
code or a distance 2 code both of which have an upper bound. Keeping track of 
the change in the parameters will give us an upper bound on the parameters of the 
original code. 

Let w = (ai,a 2 , . . . , a n \h, b 2 , ■ ■ ■ ,b n ) G F 2n . We denote by p(w) G F 2 "" 2 , the 
vector obtained by deleting the first and the n+l th coordinates of w. Thus we have 

p(w) = (o2, • • • , a n \b 2 , ...,&„) G IF 2 ™- 2 . 

Similarly, given a classical code C C F 2 ™ we denote the puncturing of a codeword 
or code in the first and n + 1 coordinates by p(C) . 

For F g -linear codes instead of considering the trace symplectic inner product we 
can consider the relatively simpler symplectic product. The symplectic product of 
u = (a\b) andv = (a'\b') in F 2 ™ is defined as (u\v) s = ((a\b)\(a'\b')) s =a'-b-a-b'. 
The symplectic dual of a code C C F 2 ™ is defined as C ±3 = {x e F 2 ™ | (x\y) s = 
0, for all y G C}. It will be seen that (u\v) t — tr q / p ((u\v) s ). 

Lemma 3.1. Let C C F 2 ™ be an ¥ q -lin ear code with (a\b) € C and (a'\b') G C^* . 
Then ((a|6)|(a / j6')>* = if and only if ((a\b)\(a'\b')) s = a ■ b' - a' ■ b = 0. It follows 
that C u = . 

Proof. If ((a\b)\(a'\b')) s — 0, then tr g / p (a' • b — a - b') ~Q. Since C is linear (aa\ab) 
is also orthogonal to (a'\b') for any aeFJ. Hence, tr q / p (aa' ■ b — aa ■ b') = 0. But 
tr is a nondegencratc function. It follows that a' • b — a ■ b' = 0. The converse is 
straightforward. The equality of C^* = C ±B follows immediately from the first part 
of the statement. □ 

As we shall be concerned with F g -lincar codes in this paper, we will focus only 
on the symplectic inner product in the rest of the paper. 

Lemma 3.2. Let C C F 2 ™ be an ¥ q -linear code. Then C has an ¥ q -linear basis of 
the form 

B = {Zl, . .. ,Zk, Zk+l,Xk+l,Zk+2, Xk+2, ■ ■ ■ , Zk+r,Xk+r} 

where (xi\xj) s = = (zi\zj) s and (xi\zj) s = 5ij. 

Proof. First we choose a basis B — {z\, . . . , Zk} for a maximal isotropic subspace 
Co of C. If Co 7^ C, then we can choose a codeword xi in C that is orthogonal to 
all of the Zk except one, say Z\ (renumbering if necessary). We can scale x\ by an 
element in F^ so that (zi|a;i} s = 1. If (Co,a;i) ^ C, then we repeat the process 
until we have a basis of the desired form. □ 

For the remainder of the section, we fix the following notation. By theorem 2.1, 
we can associate with an F g -linear [[n, k, r, d]] q subsystem code two classical F g - 
linear codes C, D C F 2 " such that D = C f~\ C ±s , |C| = q n - k + r , \D\ = g «- fe -'' an d 
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swt(Z? ±s \C) — d. By lemma 3.2, we can also assume that C is generated by 

C (zi , . . . , z s: z s +i , x s +i , . . . , z s j rT: x s + r ) , 

where s = n — k — r and the vectors Xi, z% in ¥ q n satisfy the relations (xi\xj) s = 
= (zi\zj) s and (x,|,Zj) s = <5jj. These relations on Xi,Zi imply that 

! z s , s _|_ r _|_i , a; s _|_ r _|_i, • ■ • , Zg+r+k, a; s + r +fe)i 

1 Zs > ^3+1 ) ^S+l ) • ■ • ) ^71 , X n ) . 

Lemma 3.3. An ¥ q -linear [[n,k,r > 0, d > 2]] q Clifford subsystem code exactly 
pure to 1 can be punctured to an ¥ q -linear [[n — 1, fc, r — 1, > code. 

Proof. As mentioned above, we can associate to the subsystem code two classical 
codes C, D C Fg™. Two cases arise depending on swt(D). 

a) If swt(-D) = 1, then without loss of generality we can assume that swt(zi) = 1. 
Further, z\ can be taken to be of the form (1,0,..., 0\a, 0, . . . , 0). and for i ^ 1, 
because of F 9 -linearity of the codes we can choose every Xi, Zi to be of the form 
(0, a 2 , . . . , a n \bi, &2> • • • ) b n ). Further, as Xi, Zi must satisfy the orthogonality 
relations with z\ viz., {z\\z^) s = = (zi\xi) s , for i > 1 we can choose x^, Zi to 
be of the form (0, a 2 , . . . , o„|0, 6 2 , . . . , b n ). It follows that because of the form 
of Xi and z% puncturing the first and n + 1 th coordinate will not alter these 
orthogonality relations, in particular (p(xi)\p(zi)) s ^ for s + 1 < i < n. 
Letting p{xi) = x\, p{zi) = z[ and observing that p{z\) = (0, . . . , 0|0, . . . , 0), 
we see that the code p(C) = (z' 2 , . . . , z' s , z' s+1 , x' s+1 , . . . , z' s+r , x' s+r ). Denoting 
by Dp — p(C) n p(C) ±s it is immediate that D p is generated by {z' 2l . . . , z' s } 
while Dp B = (z 2 ,...,z' s ,z' s+1 ,x' s+1 ,...,z' n ,x' n ). Hence p(C) defines an [[n - 
l,fc,r,swt(D^ s \p(C))]] q code. 

Next we show that swt(Dp s \ p{C)) > d. Let u — (a 2 , . . . , a n \b 2 , . . . ,b n ) be 
in Dp B \ p(C), then we can easily verify that (0, a 2 , . . . , a„|0, b 2 , . . . , b n ) is 
orthogonal to all z^, 1 < i < s and hence it is in D- 1 '. It cannot be in C 
as that would imply that u is in p(C). But swt(D ± ' \ C) > d. Therefore 
swt(u) > d. and p(C) defines an \[n — l,k,r,> d]] q code. By choosing C = 
{z 2 , . . . , z' s , z' s+1 , z' s+2l x' s+2 , . . . , z' s+rl x' s+r ) wc can conclude that there exists 
an [[n— , k, r— 1, d]] q code. Alternatively, apply theorem 16 in Aly et al. (2006). 

b) If swt(D) > 1, then we can assume that swt(z s +i) = 1 and form the code 
C = (zi, . . . ,z s ,z s+ i,z s+2 ,x s+2 , . . . ,z s+r ,x s+r ). It is clear that C defines 
an [[n, k, r — 1, d]] q code that is pure to 1 with swt(C" fl C" ±s ) = 1. But this 
is just the previous case, from which we can conclude that there exists an 
[[n — 1, k, r — 1, > d]] q code. 

□ 

Lemma 3.3 allows us to establish a bound for distance 2 codes which can then 
be used to prove the bound for arbitrary distances. 

Lemma 3.4. An impure ¥ q -linear [[n, k,r,d = 2]] q Clifford subsystem code satisfies 

k + r<n-2d + 2. 
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Proof. Suppose that there exists an F^-linear [[n, k,r, d = 2]] q impure subsystem 
code such that k + r > n — 2d + 2; in particular, this code must be pure to 1. By 
lemma 3.3 it can be punctured to give an [[n—l,k,r — l,> d]] 2 subsystem code. If 
this code is pure, then k + r — 1 <n — 1 — 2rf + 2 holds, contradicting our assumption 
k + r > n — 2d +2; hence, the resulting code is once again impure and pure to 1. 

Now we repeatedly apply lemma 3.3 to puncture the shortened codes until we 
get an [[n — r, k, 0, > d]] q subsystem code. But this is a stabilizer code which must 
obey the Singleton bound k < n — r — 2d +2, contradicting our initial assumption 
k + r > n — 2d +2. Therefore, we can conclude that k + r < n — 2d +2. □ 

If the codes are of distance greater than 2, then we puncture the code until 
it either has distance 2 or it is a pure code. The following result tells us how the 
parameters of the subsystem codes vary on puncturing. 

Lemma 3.5. An impure ¥ q -linear [[n, k,r,d> 3}] q Clifford subsystem code exactly 
pure to d' > 2 implies the existence of an ¥ q -linear [[n — 1, k, r + 1, > d — l]] q 
subsystem code. 

Proof. Recall that the existence of an [[n, k,r,d > 3]] q subsystem code implies the 
existence of F g -linear codes C and D such that 

with s = n — k — r, and D = C n C ±3 , see above. 

The stabilizer code defined by D satisfies k + r = n — s < n — 2d + 2, or 
equivalently s > 2d — 2; it follows that s > 2, since d > d' > 2. Without loss of 
generality, we can take z\ to be of the form (1, a-2, ■ ■ ■ , a n \b\, hi . . . , b n ) for if no such 
codeword exists in D, then (0, 0, . . . , 0| 1, 0, . . . , 0) is contained in D ±B , contradicting 
the fact that swt(_D J - 3 ) > 2. Consequently, we can choose Z2 in D to be of the form 
(0, C2, • • • , c„|l, d2, . . . , d n ), and we may further assume that b\ = in z\. The form 
of z\ and Z2 allows us to assume that any remaining generator of C is of the form 
(0,« 2 , • • • ,u n \0,v 2 , ■ • • , v n ). 

Let p be the map defined by puncturing the first and (n + l) th coordinate of a 
vector in C. Define for all i the punctured vectors x[ = p{xi) and z[ — p(zi). Then 
one easily checks that (p(xi) | p(xj)) s = = (p(zi) | p(zj)) s for all indices i and j, 
and (p(xi) | p(zj)) s — Si,j if i > s + 1 or j > 3, and that {p{z\) \ p{z2)) s = — 1- 

Let us look at the punctured code p(C), 

p(C) = (z 3 , . . . , z s , z s+ i, x s+ i, . . . , z s+r , x s+rl z 1 , z 2 ). 

Since (p(z\) | p{z2)) s = —1 we have D p = p(C) n p(C) ±s — (z' 3 , . . . , z' s ), whence 
\D p \ = \D\/q 2 . As swt(C) > 2, it follows that \p(C)\ = \C\. Thus p(C) defines an 
[[n — 1, k, r + 1, swt(Dp s \ p(C))]] q subsystem code. 

Recall that the code D is generated by s > 2 vectors; we will show next that our 
assumptions actually force s > 3. Indeed, if s = 2, then \D\ = q 2 and \D ±S | = q 2n ~ 2 . 
Under the assumption swt(L>- Ls ) > 2, it follows that ^(D^OI = \ D±3 \ = q 2n ~ 2 - 
But as p{D^) C F 2n ~ 2 this implies that piD^-*) = ¥ 2 q n - 2 . Since ¥ 2n ~ 2 has 2n - 2 
independent codewords of symplectic weight one, D ±s must have 2n— 2 independent 
codewords of symplectic weight two. However, this contradicts our assumptions on 
the minimum distance of the subsystem code: 



Article submitted to Royal Society 



8 



Klappenecker, Sarvepalli 



(a) If C is a proper subspace of D ±s , then the minimum distance d is given by 
d = swt(D ±s \ C) > 3; thus, the weight 2 vectors must all be contained in C, 
which shows that \C\ = q 2n ~ 2 = \D\, contradicting |C| < |£> ±s |- 

(b) If C = D ±s , then the minimum distance is given by d — swt(D ±s ) = 2, 
contradicting our assumption that d > 3. 

Thus, from now on, we can assume that s > 3. 

Before bounding the minimum distance of the punctured subsystem code, we are 
going to show that D^ s = p(D ±s ). Let w = (ui,u 2 , ■ ■ ■ ,u n \vi,V2, ■ ■ ■ ,v n ) be a vec- 
tor in D^- s . For 3 < i < s, the vectors Zi are of the form (0,ft2, . . . ,a„|0,&2 ; • • • ,b n ); 
thus, it follows from (w\zi) s — that (p(w)\z' i ) s = 0. Hence p(w) is in Dp", which 
implies p(D^) C D£' . We have \D£'\ = q 2n - 2 /\D p \ = q 2n /\D\ = \D^\, and we 
note that {D^l = \p(D J -")\, because swt(L»- L =) > 2; hence, = p{D^°). 

Let w' = (u2, . . . ,u n \v2, ■ ■ ■ ,v n ) be an arbitrary vector in p(D- Ls ) \ p(C). It 
follows that there exist some a, (3 in ¥ q such that w = (a, ti 2 , ■ • ■ , u n \/3, v 2 , ■ . ■ , v n ) 
is in D ±s ; it is clear that w cannot be in C, since then p(w) — w' would be in p(C); 
hence, swt(w) > d. It immediately follows that swt(D^ s \ p(C)) > d — 1. Hence 
p(C) defines an [[n — 1, k, r + 1, > d — l]] q subsystem code. □ 

Now we are ready the prove the upper bound for an arbitrary subsystem code. 
Essentially we reduce it to a pure code or distance two code by repeated puncturing 
and bound the parameters by carefully tracing the changes. 

Theorem 3.6. An ¥ q -linear [[n, k,r, d > 2]] q Clifford subsystem code satisfies 

k + r<n-2d+2. (3.1) 

Proof. The bound holds for all pure codes, see Aly et al. (2006). So assume that 
the code is impure. If d = 2, then the relation holds by lemma 3.4; so let d > 3. If 
the code is exactly pure to 1, then it can be punctured using lemma 3.3 to give an 
[[n — 1, k, r — l,d' — d]] q code, otherwise it can be punctured using lemma 3.5 to 
obtain an [[n — 1, k, r + 1, d' > d — l]] q code. If the punctured code is pure, then 
it follows that either fc + r- l<n-l-2d+2orfc + r + l<n-l-2<i' + 2< 
n — 1 — 2(<i— 1) + 2 holds; in both cases, these inequalities imply that k+r < n — 2d+2. 

If the resulting code is impure, then if it is exactly pure to 1 we puncture the code 
again using lemma 3.3, if not we puncture using lemma 3.5, until we get a pure code 
or a code with distance two. Assume that we punctured i times using lemma 3.3 and 
j times using lemma 3.5, then the resulting code is an [[n—i—j, k, r+j—i, d' > d— j]] g 
subsystem code. Since pure subsystem codes and distance 2 subsystem codes satisfy 

k + r + j - i < n - i - j - 2d' + 2 < n - i - j - 2{d - j) + 2, 

it follows that k + r<n-2d+2 holds. □ 

When the subsystem codes are over a prime alphabet, this bound holds for all 
codes over that alphabet. In the more general case where the code is not linear, 
numerical evidence indicates that it is unlikely that the additive subsystem codes 
have a different bound. We have shown that a large class of impure codes already 
satisfy this bound. We conjecture that all subsystem codes satisfy k + r < n — 2d+2. 
Next, we give an application of this upper bound. 



Article submitted to Royal Society 



Subsystem Codes 



9 



(b) Can subsystem codes improve upon MDS stabilizer codes? 

In this subsection, we compare stabilizer codes with subsystem codes. We first 
need to establish the criteria for the comparison, since subsystem codes cannot be 
universally better than stabilizer codes. For example, it is known that an [[n, k, r, d]] q 
subsystem code can be converted to an [[n, k,d]] q stabilizer code (see Aly et al. 
(2006), lemma 10 for a proof of this claim); this implies that no [[n, k, r, d]] q sub- 
system code can beat an optimal [[n, k, d']] q stabilizer code in terms of minimum 
distance, as d' > d. One of the attractive features of subsystem codes is a potential 
reduction of the number of syndrome measurements, and we use this criterion as 
the basis for our comparison. 

First, we must highlight a subtle point on the required number of syndrome bits 
for an F g -lincar [n, k, d] q code. A complete decoder, will require n — k syndrome 
bits. Complete decoders are also optimal decoders. A bounded distance decoder on 
the other hand can potentially decode with fewer syndrome bits. Bounded distance 
decoders typically decode up to [(d— 1)/2J. However, to the best of our knowledge, 
except for the lookup table decoding method, all bounded distance decoders also 
require n — k syndrome bits. As the complexity of decoding using a lookup table 
increases exponentially in n— k it is highly impractical for long lengths. We therefore 
assume that for practical purposes, that we need n — k syndrome bits. 

Similarly, for an F 9 -linear [[n, k, r, d]] q subsystem code, a complete decoder will 
require n — k — r syndrome measurements, as is shown in Appendix A. We are 
not aware of any quantum code, stabilizer or subsystem, for which there exists a 
bounded distance decoder that uses less than n — k — r syndrome measurements to 
perform bounded distance decoding. The work by Poulin (2005) prompts the follow- 
ing question: Given an optimal [[k+ 2d— 2, k, d]] q MDS stabilizer code, is it possible 
to find an [[n, k, r, d]] q subsystem code that uses fewer syndrome measurements? 

There exist numerous known examples of subsystem codes that improve upon 
nonoptimal stabilizer codes. The fact that the stabilizer code is assumed to be 
optimal makes this question interesting. The Singleton bound k + r < n — 2d + 2 
of an F 9 -linear [[n, k, r, d]] q subsystem code implies that the number n — k — r of 
syndrome measurements is bounded by n — k — r > 2d — 2; thus, for fixed minimum 
distance d, there exists a trade off between the dimension k and the difference n — r 
between length and number of gauge qudits. 

Corollary 3.7. Under complete decoding an ¥ q -linear [[n 7 k,r,d > 2]] q Clifford 
subsystem code cannot use fewer syndrome measurements than an ¥ q -linear [[k + 
2d — 2, k, d]] q stabilizer code. 

Proof. Seeking a contradiction, we assume that there exists an [[n, k, r, d]] q sub- 
system code that requires fewer syndrome measurements that the optimal \[k + 
2d — 2, k, d]] q MDS stabilizer code. In other words, the number of syndrome mea- 
surement yield the inequality k + 2d — 2 — k > n — k — r, which is equivalent to 
k + r > n — 2d +2, but this contradicts the Singleton bound. □ 

Poulin (2005) showed by exhaustive computer search that there does not exist 
an [[5, 1, r > 0, 3]]2 subsystem code. The above result confirms his computer search 
and shows further that not even allowing longer lengths and more gauge qudits can 
help in reducing the number of syndrome measurements. In fact, we conjecture that 
corollary 3.7 holds for bounded distance decoders also. 
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We wish to caution the reader that gains in error recovery cannot be quantified 
purely by the number of syndrome measurements. In practice, more complex mea- 
sures such as the simplicity of the decoding algorithm or the resulting threshold 
in fault-tolerant quantum computing are more relevant. The drawback is that the 
comparison of large classes of codes becomes unwieldy when such complex criteria 
are used. 

4. Subsystem codes on a lattice 

Bacon gave the first family of subsystem codes generalizing the ideas of Shor's 
[[9, 1, 3]]2 code (Bacon 2006). Recently, he and Casaccino gave another construction 
which generalizes this further by considering a pair of classical codes (Bacon & 
Casaccino 2006). We show that this method is a special case of theorem 2.1. Since 
this construction is not limited to binary codes and our proofs remain essentially 
the same, we will immediately discuss a generalization to nonbinary alphabets. 

Theorem 4.1. For i G {1,2}, let Ci C F™ ; be ¥ q -linear codes with the parameters 
[rii, ki, di] q . Then there exists a Clifford subsystem code with the parameters 

[[nm 2 , fak 2 , (m - ki)(n 2 - k 2 ), min{di, d 2 }]] q 

that is pure to d p — minjd^, d^}, where df denotes the minimum distance of Ci . 

Proof. Let C be the classical linear code given by C = (F^ 1 ® Ci) x (Ci <8> F£ 2 ). 
Then dimC = n\(n 2 — k 2 ) + n 2 (n\ — k\) and swt(C \ {0}) > min{di, d^}. The 
symplectic dual of C is given by 

C ±s = (Ci ® F£ 2 ) x x (F™ 1 ® ci) 1 - 
= (Ci <8> F™ 2 ) x (F™ 1 <x> C 2 ). 

We have dimC ±s = k\n 2 + n\k 2 . The code D — C D C ±s is given by 
D = ((F^ 1 ® Ci) x {Ci ® F£ 2 )) n ((Ci ® F r g 12 ) x (F; 11 ® C 2 )) 

= {(F r q L1 ® c£) n (Ci ® f; 12 )) x ({Ci ® f; 12 ) n (F^ 1 ® c 2 )) 

= (Ci®C 2 L ) x (C^Ca), 

and dimD = k\(n 2 — k 2 ) + k 2 {n\ — k\). It follows that dimC — dim!) = 2(ni — 
fci)(«2 — fo) and dimC ±a — dimD = 2k\k 2 . Using corollary 2.2, we can get a 
subsystem code with the parameters 

[[n 1 n 2 ,k 1 k 2 , (m - fci)(n 2 - k 2 ), d = swt(L>- Ls \ C)]] q 

that is pure to d p — min-fc^, d 2 }. It remains to show that d = mm{di,d 2 }. 
Since D = (Ci ® C 2 ) x {Ci ® C 2 ), we have 

£ ±s = {Ci ® C2) 1 - x (Ci ® Ci) 1 - 

= ((d ® F™ 2 ) + (F™ 1 ® C#)) x ((F™ 1 ® C 2 ) + (C^ ® F£ 2 )) . 

In the last equality, we used the fact that vectors u\ ®u 2 and v\ ®v 2 are orthogonal 
if and only if u\ _L v\ or u 2 -L v 2 . 
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For i G {1, 2}, let Gi and Hi respectively denote the generator and parity check 
matrix of the code Cj. Without loss of generality, we may assume that these matrices 
are in standard form 

H t = [ I ni - ki P t ] and d = [ -P\ I ki ] , 

where P* is the transpose of Pi. Let H f = [ I ki ]■ Using these notations, the 
generator matrices of C and D ±s can be written as 



G c 









Hi® I n 



and G D ± 



Gi ® Hi 
I ni ® H 2 







H x 



Go 



It follows that the minimum distance d is given by 



swt(D" Ls \ G) = min \ wt 



wt 



Gi ® F| 
® H 2 

Hi ® G 2 



\ < / ni (8) ^2 ) 



\< ^1 



) 



Let us compute 



wt 



\< ^1 



) 



J?f <8> G 2 

If minimum weight codeword is present in D^ s \ G, it must be expressed as linear 
combination of at least one row from [iJf (g> G 2 \ otherwise the codeword is entirely 
in G. Recall that H\ = [ I ni -k 1 Pi ] and H\ = [ /fej ]. Letting Pi = (pij), 
we can write 












. 


. 


G 2 














. 


. 





G 2 


Hit. 


3G 2 " 







. 


. 





G 2 


Hi (5 






I 712 


. 


. 




Plk^n.2 















P2lln 2 


PlkJm 











. 


-'712 


P(m-fei)l-^n 2 


P(m — ki)ki^n2 



Now observe that any row below the line in the above matrix can has a weight of 
only one in each of the last k\ blocks of size n 2 . And any linear combination of them 
involving less than c? 2 and at least one generator from the rows above must have a 
weight > d 2 . If on the other hand there are more than d 2 rows involved, then the 
first n 2 (ni — fci) columns will have a weight > d 2 . Thus in either case the weight 
of an element that involves a generator from [iff eg) G 2 ] must have a weight > d 2 . 
On the other hand, the minimum weight of the span of [H^ <g> G 2 ] is wt(G 2 ) = d 2 , 
from which we can conclude that 



wt 



H{ ® G 2 
Hi ® I n2 



) 



d 2 
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Because of the symmetry in the code we can argue that 



wt 



di 



and consequently d = min{di, g?2j, which proves the theorem. 



□ 



(a) Bacon-Shor codes 

Bacon (2006) proposed one of the first families of subsystem codes based on 
square lattices. A trivial modification using rectangular lattices instead of square 
ones gives the following codes, see also Bacon & Casaccino (2006). The relevance 
of these codes will be seen later in §5. Using the same notation as in theorem 4.1, 
let Gi — [1, . . . , l]i X j and Hi be the matrix defined as 



Hi 



1 1 
1 1 



1 1 
1 1 



-1X2 



and C, the additive code generated by the following matrix. 



G 



I ni ® H n2 





H ni <S> I n 



Observe that Gi generates an [i,l,i] q code with distance i. By theorem 4.1, G ni 
and G„ 2 will give us the following family of codes 

Corollary 4.2. There exist [[nin 2 , 1, (ni — l)(n 2 — 1), min{ni, n 2 }]] 9 Clifford sub- 
system codes. 



5. Subsystem codes and packing 

We investigate whether subsystem codes lead to better codes because of the decom- 
position of the code space. Since the early days of quantum codes, it has recognized 
that the degeneracy of quantum codes could lead to a more efficient quantum code 
and allow for a much more compact packing of the subspaces in the Hilbert space. 
But so far it has not been shown for stabilizer codes. We can derive similar bound 
for subsystem codes. Aly et al. (2006) showed the following theorem for pure sub- 
system codes. 

Theorem 5.1. A pure ((n,K,R,d)) q Clifford subsystem code satisfies 

L(d-l)/2j 

E L )(l 2 -l) j <Q n /KR. (5.1) 

It is natural to ask if impure subsystem codes also satisfy this bound. We 
show that they do not by giving an explicit counterexample. This counter example 
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comes from the codes proposed by Bacon (2006). Recall the Bacon-Shor codes are 
[[n 2 ,l,(n — l) 2 ,n]]2 subsystem codes. The [[9,1,4, 3]]2 is an interesting code. We 
can check that it satisfies the Singleton bound for subsystem codes as 



So it is an optimal code. More interestingly, substituting the parameters of the 
[[9, 1, 4, 3]]2 Bacon-Shor code in the above inequality we get 



Therefore the [[9, 1, 4, 3]]2 Bacon-Shor code beats the quantum Hamming bound for 
the pure subsystem codes proving the following result. 

Theorem 5.2. There exist impure ((n, K, R,d)) q Clifford subsystem codes that do 
not satisfy 



An obvious question is why impure codes can potentially pack more efficiently 
than the pure codes. Let us understand this by looking at the [[9, 1,4, 3]]2 code a 
little more closely. This code encodes information into a subspace, Q where dimQ = 
2 k + r = 2 5 . As it is a subsystem code Q can be decomposed as Q = A eg) B, with 
dimA = 2 fe = 2 and dim_B = 2 r = 2 4 . In a pure single error correcting code all 
single errors must take the code space into orthogonal subspaces. In an impure code 
this is not required two or more distinct errors can take the code space to the same 
orthogonal space. In the Bacon-Shor code a phase flip error on any of the first three 
qubits will take the code space to same orthogonal subspace and because of this 
we cannot distinguish between these errors. However, it is not a problem because 
we can restore the code space with respect to A even though we cannot restore B. 
Thus instead of requiring 9 orthogonal subspaces as in a pure code, we only require 
3 orthogonal subspaces to correct for any single phase flip error. Considering the bit 
flip errors and the combinations we need only 9 orthogonal subspaces. Thus with 
the original code space this means we need to pack ten 2 5 -dimensional subspaces 
in the 2™ = 2 9 dimensional ambient space, which is achievable as 10 • 2 5 < 2 9 . 

More generally, in a sense degeneracy allows distinct errors to share the same 
orthogonal subspace and thus pack more efficiently. It must be pointed out though 
that this better packing is attained at the cost of r gauge qudits compared to a 
stabilizer code. 

In fact there exists another code among the Bacon-Shor codes which also beats 
the Hamming bound for the subsystem codes. This is the [[16, 1,9,4]] 2 code. The 
family of codes given in corollary 4.2 provides us with [[12, 1, 6, 3]]2, yet another 
example of a code that beats the quantum Hamming bound like the [[9, 1,4, 3]]2 
code. We can check that 



fc + r = l + 4 = n-2rf + 2 = 9- 6 + 2. 




E 
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But note that unlike [[9, 1,4, 3]]2 this code does not meet the Singleton bound for 
pure subsystem codes as 6 + 1 < 12 — 6 + 2. Naturally we can ask if there is a 
systematic method to construct codes that beat the quantum Hamming bound. At 
the moment we do not know. It appears unlikely that there exist long codes that 
beat the quantum Hamming bound. 



6. Conclusion 

We have proved that any F ? -linear [[n, k, r, d]] q Clifford subsystem code obeys the 
Singleton bound k + r < n — 2d + 2. Furthermore, we have shown earlier that 
pure Clifford subsystem codes satisfy this bound as well. Our results provide much 
evidence for the conjecture that the Singleton bound holds for arbitrary subsystem 
codes. 

Pure Clifford subsystem codes obey the Hamming (or sphere packing) bound. 
In this paper, we have shown the amazing fact that there exist impure Clifford 
subsystem codes beating the Hamming bound. This is the first illustration of a case 
when impure codes pack more efficiently than their pure counterparts. One example 
of a code beating the Hamming bound is provided by the [[9, 1, 4, 3]]2 Bacon-Shor 
code; this remarkable example also illustrates the following noteworthy facts: 

a) The [[9, 1, 4, 3]]2 code requires 9 — 1 — 4 = 4 syndrome measurements just like 
the perfect [[5, 1,3]]2 code. 

b) Since k + r < n — 2d + 2 for all prime alphabet codes, [[9, 1, 4, 3]]2 code is also 
an optimal subsystem code. This is interesting because the underlying classical 
codes are not MDS. In MDS stabilizer codes, the underlying classical codes are 
required to be MDS codes. 

c) The Bacon-Shor code is also impure. So unlike MDS stabilizer codes which must 
be pure, MDS subsystem codes can be impure. 

d) The maximal length of a g-ary stabilizer MDS code is 2q 2 — 2, (Ketkar et al. 
2006) whereas for subsystem codes it is larger as the [[9, 1, 4, 3]]2 code indicates. 

The implication of b)-d) is that optimal subsystem codes can be derived from 
suboptimal classical codes, unlike stabilizer codes. 

We conclude with a few open questions that seem worth investigating. 

i) Do arbitrary [\n, k, r, d]] q subsystem codes also satisfy k + r < n — 2d +21 

ii) Is the Hamming bound for subsystem codes obeyed asymptotically? 

iii) What is the maximal length of MDS subsystem codes? 

The second question is motivated by the fact that binary stabilizer codes obey the 
quantum Hamming bound asymptotically, see Ashikhmin & Litsyn (1999). 
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Appendix A. Syndrome measurement for nonbinary 

Fg-linear codes 

Decoding of nonbinary quantum codes has not been studied as well as binary codes. 
Encoding of F g -linear nonbinary quantum codes was investigated in Grassl et al. 
(2003). The authors suggest that the decoder is simply the encoder running back- 
wards. While that maybe reasonable in quantum communication, it is not preferable 
in the case of quantum computation. 

Here we give a method that allows us to measure the syndrome for F 9 -linear 
nonbinary quantum codes. We also show that an F 9 -linear [[n, k, r, d]] q code requires 
n — k — r syndrome measurements. But first we need the definition of the following 
nonbinary gates, see Grassl et al. (2003). 

i) X(a) \x) = \x + a) 

ii) Z(b) \x) = uj tr ^/p^ \x),w = e^'P 

iii) M(c) \x) = \cx),cG F* 

iv) f \x) = ^ q Y Jy&q ^ (xv) \y) 

v) A\x) \y) = \x) \x + y) 

Graphically, these gates are represented below. 



— X(a) — — Z(b) — 



— F — 



i) ii) 
Consider the following circuit. 

\a) - 

\y) — iff? 



in) 



-e- 



iv) 
\a) 

\y + ag x ) 



-e- 

v) 



Alternatively, this circuit maps \a) \x) to \a) X(ag x ) \y). Observe that this circuit 
effectively applies X(ag x ) on the second qudit. Using the linearity, we can analyze 
the following circuit. 

|0) - 



\y) 



-0- 



E Q eF I") \y + a 9x) 



The above circuit maps |0) \y) to EaeF, \ a ) X{ag x ) \y). Using the fact that FX(b)F^ = 
Z(b), we can show that the following circuit maps |6) \y) to \b) Z(bg z ) \y). 



\b) 

\y) 



Ft 



91 



Z(bg z ) \y) 



If we wanted to apply a general operator X(ag x )Z(ag z ) to the second qudit condi- 
tioned on the first one, then we can combine the previous circuits as follows. 



\a) 

\y) 



Ft 



B— 9 



®— 9 



\a) 

X(ag x )Z(ag z ) \y) 



The above implementation is not optimal in terms of gates, but it will suffice for our 
purposes. Consider an [[n, k, r, d]] q code. Let E be an error in £ . If E is detectable, 
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then E does not commute with some element (s) in the stabilizer of the code. Let 

9 = {9x\gz) = {0, . . . ,0, aj , . . . ,a n \0, . . . ,0,bj, . . . ,b n ) G Ff , 

where (oj, 6j) ^ (0, 0), be a generator of the stabilizer. Then for all detectable errors 
that do not commute with a multiple of g, the following circuit gives a nonzero value 
on measurement. 



|0> 



F 



-j Ft - 6. 1 



M -L Ft - &" 1 



bj — F 



-©- 



b n — F 



Ft 



-e- 



Note that whenever (oj, 6j) = (0, 0), then we leave that qudit alone. Similarly if a, 
or &i are zero, then we do not implement the corresponding portion. Let the input 
to the above circuit be E \ where is an encoded state. It can be easily verified 
that the above circuit maps the state |0) E \tp) to 

^ \a) X(ag x )Z(ag z )E . 

qGF, 

Let X{g x )Z{g z )E = cj tT "/p^EX(g x )Z(g z ), where X(g x )Z(g z ) is corresponding ma- 
trix representation of g. Then we have X(ag x )Z(ag z )E — uj tI "/p^ at " > EX(g x )Z(g z ), 
by lemma 5 in (Ketkar et al. 2006). Thus we can write 

J2\a)X(ag x )Z(ag z )E\i>) = £ \a) co^^ EX(ag x )Z(ag z ) |V>) , 

u£F q a£F q 



\a)u tT «'pW F| 



K a£F„ 



where we have made use of the fact that X(ag x )Z(ag z ) 
is in the stabilizer. The final state is given by 



as X(ag x )Z(ag z ) 



F^\a)X{ a g x )Z{ a g z )E\^) 

a£F q 



J2 \a) u tT «>" {at) E |V) , 

a£F q 

^ ^ uj-^/pW) \/3) uj tT */»( at) E 

a£F q f3<EF q 
0£F q a£F q 

53 1/3) 53 ^ tIq/p{at ^ al3) E\^), 

fi£F q a£F q 

\t)EW), 
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where the last equality follows from the property of the characters of ¥ q . Next 
we observe that the error aE, where a £ ¥ q gives \at) on measurement. Strictly 
speaking we refer to the preimage of aE in £ . Hence the syndrome qudit can take 
q different values. Since every detectable error docs not commute with some IF g - 
multiple of a stabilizer generator, we have the following lemma on the necessary 
and sufficient number of syndrome measurements. 

Lemma 6.1. Given an ¥ q -linear [[n, k, r, d]] q Clifford subsystem code, n — k — r 
syndrome measurements are required for decoding it completely. 

Proof. Let g be a generator of the stabilizer of the subsystem code. By theorem 2.1 
and lemma 3.2, for every generator g there exists at least one detectable error that 
does not commute with g but commutes with all the other generators. This error 
can be detected only by measuring g. Thus we need to measure all the generators 
of the stabilizer, equivalcntly n — k — r syndrome measurements must be performed. 

Every correctable error takes the code space into a q' c+r -dimcnsional orthogonal 
subspace in the (/"-dimensional ambient space, see §2. Each of these errors will give a 
distinct syndrome. This implies that we can have g n ^ fe ^ r distinct syndromes. Since 
each syndrome measurement can have q possible outcomes and there are n—k — r 
generators, these measurements are sufficient for performing error correction. □ 

This parallels the classical case where an [n, k, d] q code requires n—k syndrome 
bits. A subtle caveat must be issued to the reader. If we choose to perform bounded 
distance decoding, then it maybe possible that the set of correctable errors can 
be distinguished by a smaller number of syndrome measurements. But even in the 
case of (classical) bounded distance decoding it is often the case that we need to 
measure all the syndrome bits. 
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